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ABSTRACT 

Intra Cluster Media (ICMs) located at galaxy clusters is in the state of hot, 
tenuous, magnetized, and highly ionized X-ray emitting plasmas. This overall 
collisionless, viscous, and conductive magnetohydrodynamic (MHD) turbulence 
in ICM is simulated using hyper and physical magnetic diffusivity. The results 
show that fluctuating random plasma motion amplihes the magnetic held, which 
cascades toward the diffusivity scale passing through the viscous scale. The ki¬ 
netic eddies in the subviscous scale are driven and constrained by the magnetic 
tension which hnally gets balanced with the highly damping effect of the kinetic 
eddies. However, the saturated kinetic energy spectrum is deeper than that of 
the incompressible or compressible hydrodynamics huid. To explain this unusual 
held prohle we set up two simultaneous diherential equations for the kinetic and 
magnetic energy spectrum using an Eddy Damped Quasi Normal Markovian- 
ized (EDQNM) approximation. The analytic solution tells us that the magnetic 
energy in addition to the viscous damping ehect constrains the plasma motion 
leading to the power spectra: kinetic energy spectrum Ey ~ k~^ and correspond¬ 
ing representative magnetic energy spectrum ~ Also the comparison 

of simulation results with diherent resolutions and magnetic dihusivities implies 
the role of small scale magnetic energy in dynamo. 

Subject headings: galaxies: clusters: intracluster medium, magnetic helds 


1. Introduction 

ICM located at the center of galaxy cluster is composed of fully ionized hot plasmas 
(T ~ 10®K). The gas includes most of the cluster baryons (> 85%) and heavy elements, but 
overall density ‘n’ is very low (n < 10^ cm“^). As a result ICM has very small dihusivity 
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‘?7’ (~ while viscosity V’ (~ l/'n.) is very large f Schekochihin et al.l 120051) . Dynamo 

effect from the turbulent flow motions in ICM exceeds the dissipation of such high viscous 
plasmas causing the growth of seed magnetic helds, which react back the plasma motion 
through magnetic tension (B • VB). 


The amplified magnetic field constrains s ome intrinsic properties in I CM. Non-magnetized 
thermal Spitzer conductivity is like (iNaravan fc Medvedev! l200lh : 
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where Ti = fcT/10 Kev, Ae is the mean free path of an electron and tcoui is the time between 
coulomb collisions. If magnetic field is injected, k±, conductivity perpendicular to the mag¬ 
netic field is reduced to ~ {Pe/ i^sp) so that the conductivity becomes anisotropic 
with one third of the thermal Spitzer conductivity: k = k± + k\\ —)■ k\\. This anisotropic con¬ 
ductivity brings about temperature distribution that depends on the direction and strength 
of magnetic held. Also as the conservation of magnetic moment ps = u\/B and kinetic 
energy (~ -|- u\) implies, magnetic held makes pressure tensor Pij = mmpUj anisotropic 

along the held. T he ratio of pressure aniso tropy ‘A’ is implicitly related to the viscosity and 
stability of ICM ( Schekochihin et ah 2010 ): 
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Here, kinematic viscosity z/ ~ u^f^/uu, ion-ion collision frequency uu = dvrne^ln A 
If A is smaller than —‘Ijfi or larger than 1//3, hrehose or mirror instability occurs in the 
range between ion Larmor rad ius (pi ~ lO'^ — 10® km) and the mean free path (Amfp~ 10^® 
km) ( Schekochihin et ah 2008h . Since ICM has high ‘/3’, stability range is very narrow, i.e., 
practically unstable. Instability is thought to redistribute the plasrna mo tions faster than 
coulomb collision does, but its role is not fully understood yet (iJonesI 120081) . 


By now we may begin to wonder if a typical MHD theory is applicable to the weakly colli- 
sional ICM plasma. Since collision between particles transfers momentum to make the system 
isotropic, most phenomenological or analytic methods assume the collisional MHD system 
system. However ^ since the tenuou s plasmas cannot expect sufficient collisions. Chew, Gold- 
berger, and Low (jChew et al.lll956f) developed an adiabatic MHD equation with the consider¬ 
ation of anisotropic pressure tensor (CGL-MHD model). Recently Santos -Liman insisted the 


validness o f the typical MHD theory with these additional constraints (jSantos-Lima et ah 


201ll. 120141) . Their model, a kinetic MHD based on CGL-closure with the limit of anisotropy. 
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shows the growth rate of magnetic energy is similar to that of a typical MHD theory in the 
early time regime, but the saturation is much smaller. In fact, the effect of anisotropy due 
to the magnetic field in the plasma does not appear in the kinematic regime while magnetic 
back reaction is negligibly small. 


However, regardless of the anisotropic limit, the fluids driven by a random and isotropic 
external force eventually become isotropic followed by the strongly frozen magnetic fields if 
P^M —>■ C )0 (77 —>■ 0). Even for a unit magnetic prandtl number {Ptm = vjri = 1), where 
magnetic fields are not so strongly frozen, the large scale fields driven by the random isotropic 
force te nd to remain ind ependent of direction if a background magnetic held ‘be^d is not very 
strong fjCho et al.ll2009h . Moreover the growing anisotropy in smaller scales cannot decisively 
affect the energy spectrum due to the small eddy turnover time. As long as the isotropic 
force continuously drives the system, especially large scales, the theoretical huid model that 
assumes an isotropic system is valid except the case of microscale instability or very strong 
‘bext’. In addition there is an interesting report that the magneti c energy spectrum around 
the core of cluster is possibl y Kolmogorov’s spect rum ~ k~'°P f KolmogorovI Il94ll) which 
assumes an isotropic system (IVogt fc EnBlin 11200311 . Now we can focus our interest on the 
isotropic properties in ICM of large magnetic prandtl number: energy spectra. 


Dynamo in magnetized plasma with large Ptm-, which is not rarely observed in space, occurs 
easily and has its peculiar properties. As the excited magnetic helds (energy) are instilled into 
the damped kinetic eddies, the viscous damping scale ky ~ I// 1 / is extended toward resistivity 
scale krj ~ l/lrj {ky <C kr,). Then, these two coupled velocity and magnetic held generate spe- 
cihc power spectra which are diherent from the typical spectrum of Kolmogorov’s incompress¬ 
ible huid or Burger’s c ompressible huid. There have been works on the small scale dynamo 


Yousef et ah 


for large Pr 


plasma (jSchekochihin et ahl l 2002 l: ICho et al.l 120031: ISchekochihin et al.ll2004J: 


M Pla bllla noCll6KOClllIllIl tjl Hi. I IZUUZL IL^ilO tjl ai.l IZUUOu lOCIltj 

booTh n addition to the several analytic methods JBatche 


or 


1950 


Kazantsev 


QGsl: iKulsrud and Anderson 1992 : Schober et ah 2012 : Bovino et al 2013h. To solve the 


magnetic induction equation MHD theories based on Kazantsev’s work flKazantsevi 1196 
assume the second order velocity held correlation (v ■ v) ~ {6ij — ^)T/v(r) - 1 - 
However to explain the formation of this second order velocity held correlation, i.e., ki¬ 
netic energy spectrum, we need to solve the coupled momentum and magnetic induction 
equation simultaneously instead of making an assumption of the kinetic energy spectrum in 
advance. Here we use an Eddy Damped Quasi Norrnal Ma rkovianization method EDQNM, 
(iKraichnan fc NagarajanI 119671: iMcCombI Il99fll: IParkI 120131 ) and dimensional approach in a 
limited way. We simplihed the resultant simultaneous diherential equations and found out 
the solutions for kinetic and magnetic energy spectrum. 
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In chapter 2, we briefly introduce simulation tool and analytic method used in this paper. 
Simulational and analytic results are introduced in chapter 3. And in the flnal chapter 
we discuss about the results, their physical meanings, and future topics. Detailed analytic 
calculations are discussed in the appendix. 


2. Numerical and analytic method 


We have solved the incompressible MHD equations using a pseudo-spectral code with a 
periodic box of size ‘(27r)^’: 


3 \ 
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where ‘f’ is a random mechanical force driving a system at fc ~ 2 — 3 in fourier space, 
and ‘bext’ is a weak background (guide) magnetic flelcii] {pext = 0.0001) covering the whole 
magnetic eddy scales. Here, ‘B’, magnetic held divided by ‘(dvrp)^/^’, has the unit of Alfven 
velocity, and ‘v’ is ‘rms velocity’, and ‘t’ has the unit of large scale eddy turnover time ‘L/n’. 
For example, if ‘L’ of a cluster is ~ 400 kpc and ‘n’ is ~ 400 km/s, then ‘L/n’ is ~ 10® year. 
The time ‘F has the unit of ‘10®’ year. The system realizes the state of ICM of the high 
viscous plasma state with ideally frozen magnetic fields. 


To compare the effect of hyper diffusivity (incom pressible fluid) and physical diffusivity 
(compressible fluid), we also used PENCIL CODE (IBrandenburg II 2 OOII) with message pass¬ 
ing interface(MPI) in a periodic box of spatial volume (27r)^ with mesh size 288^. The basic 
equations solved in the code are. 

Dp 
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Du 

dA 

dt 
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—— = —c^Vlnp-h 
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( 6 ) 

(7) 


^‘hext’ is not indispensable for the growth of magnetic fields if there is a substitutin g seed magnetic f ield. 


The turbulence by cosmological shocks can amplify the weak seed field of any origin (IRvu et al 


was p ointed out that subsonic turbulence could develop with a very wea k seed magnetic held (IRyu et al 


20^). Mor eover turb ulence can amplify a localized seed magnetic field (jBrandenburg 1200 1! ICho fc Yoo 


2008). It 


2012 ). Cho ( Choi 2014 ) investigated the origin of seed magnetic fields and insisted that the origin of the seed 


field should be more like the localized seed magnetic fields ejected from the astrophysical bodies 
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p: density; u: velocity; B; magnetic field; A: vector potential; J: current density; D j Dti= 
d/dt + M ■ V)\ advective derivative; p: magnetic diffusivity(=c^/47r(j, a: conductivity); z/: 
kinematic viscosity(=/r/p, p: viscosity); c^,: sound speed. Velocity is expressed in units of 
Cs, and magnetic fields in units of {pQ ^qY/‘^Cs{[B] = ^/po\J^^[v\ from Em ~ B‘^/pq and 
Ekin ~ Pov"^)- Po is magnetic permeability and po is the initial density. Note that po ~ p in 
the weakly compressible simulations. These constants Cg, po, and po are set to be ‘1’. f(a:, t) is 
represented by N fo(t) exp [ikj(t) ■x + i0(t)](A: normalization factor, fo: forcing magnitude, 
kj(f): forcing wave numbe]|§. The variables are also independent of a unit system. However, 
instead of hext, PENCIL CODE gives a system initial seed magnetic held in small scales. 
his seed held di sappears in a few simulation time steps. PENCIL CODE and the other code 


fICho et ahl 120031) are not the same in various ways, but we will see they produce practically 


the same energy spectra. 


Results 


3.1. Simulation results 


Figjl(a) shows the normalized kinetic energy Eyit) (black dashed line) and magnetic 


energy ^'^(t) (red solid line) of incompressible MHD huid {v = 0.015, p —>■ 0, Ptm —t oo, 
f3 = SttP/B"^ —)■ oo) with resolution of 256^ (thinner line) and 512^ (thicker line). Reynolds 
number Re of both cases are ~ 42, and their magnetic Reynolds number ReM are actually 
inhnity. When the random isotropic forcing begins to drive the system, Ey{t) quickly grows, 
becomes saturated, and keeps the status quo until EMif) begins to arise at f ~ 15 — 20. As 
pMit) grows, the energy transfer from Eyit) to pMit) gets accelerated. For this nonlocal en¬ 
ergy transfer, the gap between kinetic and magnetic energy is an important factor. However, 
as B ■ Vv in Eq.dH) implies, the geometrical constraint between v and B also plays a role of 
determinant. If magnetic held is parallel to the gradient of velocity held, kinetic energy is 
transferred to magnetic eddies. If not, magnetic held frozen to the plasma huid moves freely 
along the huid or gets annihilated. Around the onset position at f ~ 15 — 20, structural 
change between Ey and Em seems to get accelerated and completed at the saturation. Also 
the plot shows the onset position of EMif) is proportional to the resolution. Higher resolu¬ 
tion without resistivity makes more space where the forward cascaded EM{t) can stay not 
being dissipated much. This leads to the imbalanced energy distribution like the relatively 
smaller amount of magnetic energy in large scales and more amount of magnetic energy in 


^Pencil code selects one of 350 vectors in k/ vector set at each time step. 
\kf \ is - 1 - 2 


fo is 0.08 and injection scale 
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small scales. This makes the gap between Ey and Em in large scales, especially aronnd 
the injection scale, increase so that the nonlocal energy transfer become accelerated. The 
saturated magnetic energy Em, sat grows with the increase of resolution, but the saturated 
kinetic energy Ev,sat decreases. This is the special feature of 77 ~ 0 with a fixed high viscosity. 


Fig jl(b)| includes the normalized Ev{t) and EM{t) of compressible fluids. Their properties 
are as follows; for thicker line, PrM=75, v = 0.015, rj = 2x 10““^, Re ~ 170, Rcm ~ 1-3 x 10^, 
(3 ~ 83; for the thinner line, PrM=7500, u = 0.015, rj = 2x 10“®, Re ~ 120, ReM ~ 8.8x 10^, 
(3 ~ 430. The resolutions in both cases are 288^. Since the average of Mach number is at 
most 0.18, the effect of compressibility is not much. We can infer more Ey can be trans¬ 
ferred to Em with higher Ptm (lower rj). However, the plot shows the actual Ey^sat & Em, sat 
are inversely proportional to Ptm- Moreover Ey^^at is slightly larger than Em, sat, and this 


tendency is opposite to that of hyper diffusivity (Fig, 1(a)). If rj is small, magnetic energy 
can migrate into the smaller scales. But since the effect of dissipation grows with the wave 
number more Em is dissipated to increase energy gap between Ey and Em- This boosts 
the nonlocal energy transfer bringing about the accelerated dissipation of energy ( k^E{k)) 

in small scales. So the saturated energy Ey^^at, Em, sat are smaller than those of lower Ptm- 
However, it is not easy to conclude whether they converge on some lower limit as the physical 
diffusivity 77 —)■ 0 with the limited resolution and Ptm at this moment. 


Fig 2(a) shows the evolving energy spectrum Ey{k) and EM{k) of the incompressible fluid 
in the early time regime f ~ 0 — 10.5 (from bottom to top). And the plot in Fig 2(b) is a de 
facto the same one in the range of f = 24.9 — 32.6. When the kinetic energy begins to drive 
the system at A; = 2.5, Ey (black dashed line, A ~ 0) grows prior to Em (red solid line). 
As the advection term v ■ Vv indicates, Pm is not necessary for the local energy transfer in 
kinetic eddies. However, without Ey{k) magnetic eddies cannot receive energy from kinetic 
eddies nor transfer its energy to the neighboring eddies. Only when magnetic fields run into 
the plasma fluids with the nontrivial Ey of which gradient is in the direction of magnetic 
fields, dynamo process occurs leading to the growth of Em- Thus, in the very early time 
regime the evolution of magnetic field looks subsidiary. But around t rsj . 4 — 2.1 EM{k) in 
small scales begins to surpass Ey{k), and gets past Ey{k) which suffers from the viscous 
dissipation. So most of Ey is located near the injection scale where the viscous damping 
effect is not so much, but Em cascades forward and stays in the smaller scales. Kinetic 
eddies in k <~ 15 lose energy (or magnetic eddies receive energy through B • Vv) whereas 
smaller scale kinetic eddies {k >~ 15) receive energy through magnetic tension B ■ VB. As 
the fluid motion can amplify the magnetic field through dynamo process, also magnetic fields 
can increase the kinetic energy through Lorentz force (J x B). The magnetic fields press the 







7 


fluid through magnetic pressure (—V-B^/2) and stretch the fluid through magnetic tension. 
We will see that the unusual Ev{k) spectrum k~^ is due to not only the viscous damping 
but also the interaction with Em- 


Fig j3(a)| |3(b)| include the saturated energy spectra of incompressible fluids with hyper diffu- 
sivity. Also the saturated energy spectra of compressible fluids with physical diffusivity are 
shown in Fig,3(c), 3(d) These plots show Ey^satik) eventually converges to ~ k~^ if Pr 


M 


is not too small regardless of the different evolving prohles of Ey{t) and EM(t) in Figjl(a)| 


1(b) So the comparison of these plots will give us some clues to the formation of Ey{k) in 


small scale range. A quick look shows Ey{k) of 512^ has clear spectrum of k~^ compared 
with the kinetic energy spectrum of 256^. However, these two simulation sets have the same 
conditions: weak = 0.0001, isotropic random driving force f, injection scale kf ~ 2.5, 
viscosity u = 0.015, and negligible magnetic diffusivity rj. Also the saturated energy levels 
shown in Fig jl(a)| are not much different. Just the distribution of Em of 512^ in small scale 
is flatter than that of 256^. This implies Em in smaller scales may be a determinant of Ey 
spectrum. 


On the other hand, for the compressible system in Fig 3(c) and Fig,3(d), Ey{k) of Ptm = 
7500 has clearer and longer spectrum of k~^ despite its smaller saturated energy level: 
E{t)sat.,PrM=75 > ^(^)^at.,PrM=7500 (FiglT(^ . This means Ey^sat.{k) is not so much influ¬ 


enced by the magnitude of EM{k). However, careful look shows the peak of Em in Fig,3(d) 


is located at smaller scale regime than that of Fig, 3(c) Then kinetic eddies in wider range 
can interact with magnetic energy EM{k) whose power spectrum is less steep. In other 
words. Fig,3(b) and Fig,3(d) commonly show the slope of Em{E) is not so slanted as that of 


Fig 3(a) and Fig 3(c) 


So if we assume a critical representative spectrum EM^k) ~ k^ in the subviscous regime, 
we can infer that steeper spectrum than k"^ cannot provide the kinetic eddies with enough 
energy for Ey ~ k~^. However the exact value of ‘m’ is not known yet. Just we can guess 
this index should be negative. The scaling f actor ‘k~^\ so called an invariant scaling factor 
flRuzmaikin et al.lll982l: iKleeorin et al.lll99611 . is drawn together in next hgure for reference. 
But we will not discuss this concept further in this paper. 


Fig,4(a), 4(b) include the compensated k^Ey and k^^'^EM- Ey{k) of high Ptm has clear 
spectrum of k~^ in the subviscous scale. But for the spectrum of Em in Fig,4(b), it is 
ambiguous to pinpoint any scaling factor. Instead, we choose i.e., m = —1/2, an 


































approximate median val ue of the previous results fjCho et al.l l2003l: iLazarian et al.l I2004J: 


Sehekoehihin et al.ll2004l) . We will see analytic method derives quite exact Ev{k) with this 
scaling factor A reference line for the scaling invariant factor is drawn together. 


Fig, 5 (a) shows Em/Ev for fig, 3(a) 3(d) The ratio of hyper diffusivity is larger than that 


of physical diffusivity, which is consistent. As mentioned, with the negligible diffusivity the 
ratio proportionally depends on the resolution because of the larger Em- In contrast, for 
the case of physical diffusivity (black line), the dissipation effect growing with wave number 
(~ k'^) dissipates Em more efficiently until the energy states are balanced to be in the state 
of equipartition. 


3.2. Analytic methods & results 


3.2.1. Eddy Damped Quasi Normal Markovianization 


We start from the momentum (EqJ2D and magnetic induction equation (EqH]). Taking di¬ 
vergence of the momentum equation we can replace pressure by convection and magnetic 
tension (Leslie and Leith 1975, Yoshizawa 2011 )g: 


dvi{k,t) 

dt 


dBi{k,t) 


Migm{k)[vg{p,t)v^{r,t) - Bg{p,t)B^{r,t)] uV‘^Vi{k,t), 

p-\-r=k 

( 8 ) 

^Zmi.k)vq{p,t)Bm{r,t), (9) 

p-\-r=k 


with the definition of algebraic multipliers 


Miqm{k^ — 
Kymik) = 


2 {km^iq T kq5i, 

kqdiryi^. 


2 kq kj-^ 
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( 10 ) 


^‘v{k,ty and ^B{k,ty depend on time ‘t’ and wavenumber ‘/c’, but ‘t’ will be omitted for simplicity. 
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Then we can get the evolving second order correlation equations of {vi{k)vi{—k)) and {Bi{k)Bi{—k))\ 

A1 A2 

d ^ ^ -v 

— + 2iyk^){vi{k)vi{-k)) = ^ Mig,n(k)l(vg(p)vm(r)vi(-k)} - (vg(-p)vm(-r)vi(k)} 

p+r=k 


A3 


A4 


- (Bg(p)B^(r)vi(-k)} + (Bg(-p)Bm(-r)vi(+k)}], (11) 

B1 B2 

p+r=fc 

The third order correlation, ^Al,A2,...,B2\ is called a transport function. These triple 
correlations play a role of transfer and dissipation of energy to decide the field profiles of 
the system. If the field is helical, ‘a’ coefficient ( (j • b) - (v • uj)) for the inverse cascade 
of magnetic energy can be derived. Also other terms for the forward ca. sc ade of energies 
derived ( Krause fc Radler 1980 : Moffatt 19781: Park fc Blackman 2012a b; Pouquet et ah 


are 


19701 ). However the helical field is not assumed in our system, which means ‘a’ coefficient or 
the inverse cascade of magnetic energy is excluded. We use a quasi normalization approxi¬ 
mation, sort of an iterative method, to find the transport function (appendix). 


The formal representations of Ey{k) and Em{E) are as follows: 

= +\j dpdrei;-;(t)A(i-2y‘z‘-xyz)Evi.J>)Ev(r) 

- j dpdrei‘;,‘;{t)'^{xy + z^)Ev{r)Ev{k) - 2veEv{k) 

1 r k^ 

+- / dpdrQi;:it) — il-2y^z^-xyz)EMip)EMir) 

/ J pr 

+ / dpdre‘i^’l^it)y(y^z - z)EM(r)Evik), 


(13) 


and 


dEuik) 

dt 


= - dpdrQ'^^At)—z{l - x^)EM{r)EM{k) 


r 

„2 


dp dr QTrpi't) — {y + xz)Ev{p)EM{k) 
dp dr + xyz)Ev{p)EM{r) 


( 14 ) 
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The integral variables ‘p’, ‘r’, i.e., wavenumbers, are constrained by the relation of p + r = k. 
‘x’, ‘j/’, and ‘ 2 ;’ are cosines of the angles formed by three vectors ‘k’, ‘p’, and ‘r’. (Fig,5(a)). 
Algebraically ‘fc’, ‘p’, and ‘r’ should satisfy a condition like (Leslie and Leith 1975): 


\k — r\ < p < k + r. 


(15) 


To derive analytically solvable equations from Eq. ffTS]) . ffTTD . we need to simplify these two 
equations considering the interaction between ‘fc’ and its close wave number. So, we take 
account of only two cases: large p (A: ~ p S> r) and large r (A: ~ r S> p). In principle 
‘A:/2 ~ p ~ r’ should also be included. But since the interaction between the close wave 
vectors, local energy transfer, is dominant in the nonhelical small scale dynamo, the latter 
case can be ignored. 


With the assumption of Ev{k) ~ k"" and EM{k) ~ A:"*, we simplify the hrst term in Eq. ffTTl) 
like below: 

(i) Large p (small r, i.e., k^p^r,p^0) 

Eddy damping function Q2^p{k,t) is approximately 


QTr;ik.t) = 


_ ^-Wp'^+v{^^+r^)+IJ-kpr]t 


jyp2 p2^ 


ri=0 




{t —)■ 00 ). 


(16) 


Also as Fig,5(a), 6(a) show, we can use the relations x ~ p ~ 0, 2 ; ~ 1, and r = k — p. Then, 

rk I p2 


-I dp—^il-x^){k-prEMik) - k^EMik). 

up^ r V 


(17) 


(ii) Small p (large r, i.e., A:~rS>p, x~2:~0, p~l) 

In this case only the triad relaxation time is left. We assume that it is independent of 
time, so we can write pkpr like 

,t) - — -. (18) 

up -|- p^kpr f^kpr 



Then, 






{k 


2^m-l 


- 2A:r™ + r'”+^)x(l 


x^)EM{k) ~ 0. 


(19) 


The results, Eq. (ll7l) . (IT^ represent the first term in ‘77’ in Eq.(|2T]). The other terms can be 
found in a similar way. 
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Then, the coupled equations of Ev{k) and Euik) are 


dEv{k) 

dt 


= - {a,k^ + 2z/fc2) Evik) + {hk"^ - b2k^) 


B 


( 20 ) 


= (cifc” + C 2 *:"’+") Ev(k) - (dik'^ + 4^+") EM(k). (21) 

c'r ^^^ '<-^' 

C D 


The coefficients ‘oj’, 'bi\ 'ci\ and 'di are independent of ‘/c’, and assumed to be independent 
of time for simplicity. The matrix form of these simultaneous differential equations is simply 


■ E'v(k) - 


■ -A 

B 

' Ev{k) - 

. E'Mik) . 


C 

-D 

EmI^) 


M 


( 22 ) 


Eq. fl22l) can be solved diagonalizing the matrix ‘A^’. For this, the bases Ey{k) and Ej^ik) 
need to be transformed using a transition matrix ‘P’ which is composed of eigenvectors of 
‘Af’: 


Then, 


Ev{k) 

Ejylik) 


V 


v{k) - 

M{k) \ ’ 


V 


B B 
A + Ai A + A 2 


' V'{k) ■ 

= P~^MP 

- v{k) - 


■ Ai O' 

- v{k) - 

_ M'{k) _ 


. M{k) _ 


0 A 2 . 

_ M{k) _ 


- v{k) - 



_ M{k) _ 


_ Mo(A;)e^2i _ 


(23) 


(24) 


Ai and A 2 are eigenvalues: 

-(A + D) + ^(A + Df - AAD + ABC 

- - 2 - 

-(A + D)- ^(A + Df- - dAD + 4BC 


(25) 

(26) 


We choose a leading term in each elements with the consideration of their coefficients k ^ p 
or fc ~ r. Then the coefficients are like 


A ~ 2 i/F, B ~ C ~ C 2 /t”+^ D ~ d2k'‘*'\ 


(27) 
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Herfl since AD - BC ~ 

Ai ~ 0, A 2 ~ -2iyk^ as ‘/c’ increases. 


0 as /c —)■ 00 , the eigenvalues converge to 


Evik) and EM^k) are expressed like 


■ Ey{k) ■ 


s 

1 - 

g 

Co(A:)e~°-* 

_ Eyiik) _ 


_ 2uk‘^ 

0 

_ Mo(A:)e-2^^'* _ 


(28) 


If ‘Vo’ and ‘Mq’ are replaced by Evo{k) ~ k^ and EMoik) ~ k"^, we find the saturated 
solutions: 


EM^k) = 2i/k‘^Vo{k) ~ 2h'k 


2-E'mo(^) 


2pkP‘ 




Ey{k) = hD 


. EMojk) 

2vk? 




EMoik) ^ Evoik) 


2vk‘^ 


biD- 


^-2ukH _ ^ 


2m—2 


(29) 

(30) 


For complete energy spectra, ‘m’ is required. But it is difficult to pinpoi nt a representative 
magne tic power spectrum because Ej^ik) is a continuously changing curve. ISchekochihin et al 
(120041 ) found ‘m = O’, peak of EM{k), but we do not think a peak can be a representative 
power spectrum that drops continuously with the wavenumber ‘/c’. On the other hand, 
Lazarian et ahl (120041) found ‘m = —1’ using a simulation with a strong background mag¬ 
netic held hext, a hlling factor, and balance relation B ■ VB ~ i/V^v. So we infer the 
index ‘m = —1/2’ for the magnetic scaling factor in a system under the inhuence of a weak 
background magnetic held. Then, from Eq. ddOp we get Ey{k) ^k~^\ which matches the sim¬ 
ulation results well. Also this makes ^AD — BC ~ O’, i.e., ‘Ai ~ 0’ and ‘A 2 < 0’ so that 
the energy spectra become independent of time when they are saturated. If ‘m = —1’ is 

and EM^k) ~ k~^\ coincident with the results of (ICho et ahl 12003 


chosen, 'Ey(k' 

~ k 

Lazarian et al. 

2004) 


20041) . A simple relation Ef^/Ey = can be derived in high Ptm- 


4. Discussion 

The analytic and simulation job in this paper are to realize the high Ptm plasma state 
in ICM. Magnetic held ahects the huid motion through Lorentz force gv x B to cause the 
rotational motion of ionized particles around the magnetic held. The ehect of magnetic held, 
leading to the anisotropic system, competes with that of collision which transfers momen¬ 
tum to make the system isotropic. In fact the weakly collisional ICM plasma has few proper 


"‘The dimensional analysis of Eg. (12(11) implies 2m ^ v + 2 when the system is saturated. 
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ways to constrict the anisotropic tendency. However, if backgronnd magnetic field is not 
too strong, a system driven by a random isotropic force in large scales eventually becomes 
isotropic overall although small scale eddies still tend to be anisotropic under the influence 
of magnetic held. 


Simulation of high Ptm tells us some important features. The viscous scale is extended 
toward the much smaller diffusivity scale and nontrivial Ey in this extended scale is a 
prerequisite to the growth of Em- For the local and nonlocal energy transfer, besides energy 
gap in kinetic and magnetic eddies, additional specihc geometrical relation between v and 
B is required. For example smaller EM{k) than Ey{k) in larger scales boosts the nonlocal 
kinetic energy transfer from kinetic to magnetic eddies; and, smaller Ey{k) than EM{k) in 
small scale helps the energy transfer from magnetic to kinetic eddies. However, the energy 
transfer in magnetic eddies is possible only when B • VV (nonlocal transfer) or —V ■ VB (lo¬ 
cal transfer) is nontrivial. As a result of all these effects with the viscous damping, Ey ~ k~^ 
and smoothly changing Em are hnally saturated in the subviscous regime. Analytic analysis 
shows the viscous effect ~ coupled with E]^ induces this unusual spectrum. lEq E^ [5U]) . 


Finally, as mentioned the anisotropic features of small scale cannot affect the large scale 
driven by the isotropic force. However, if there is an instability due to the anisotropic pressure 
in microscale, plasma distribution in the whole system may change. We have not discussed 
the influence of microscale instability due to the anisotropic pressure including viscosity and 
conductivity on the MHD system in this paper, but we will leave these important topics for 
the future research on ICM. 
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A. Appendix 


For ‘Al’ in Eq. ffTT]) . we differentiate this triple correlation term over time to use Eq. 
Then, we have 

+ / + r^)] {vg{p)vm{r)vi{-k))= 

{[{^ + iyk^)vi{-k)]vg{p)vm{r)) + {vi{-k)[{^ + i^p^)vg{p)]vm{r)) 


'dt 


d 


+ {vii-k)vg{p)[{— + iyr'^)vjnir)]) = (ww) + (vvBB)... 


If we see the hrst term, for example. 


d 


(Al) 


([(^ + i'k^)vi{-k)\vj{p)vm{r))5p+r,k = ^ [Mins{-k){Vn{j)Vs{l)Vg{p)Vm{r)) 


-Mins{-k){Bn{j)Bs(l)Vq{p)Vra{r))]5j+l-k- (A2) 


the differentiation generates the fourth order correlation. Another differentiation just induces 
the fifth order correlation. So we need an assumption to close this equation. It is known 
that statistically turbulent quantities follow a normal distribution. And the fourth-order 
term {uuuu) can be decomposed into the combination of se c ond-order corr elation terms like 
below: quasi-normal approximation fjProndman et al.lll954 lTat.snmilll957l ): 


(M1M2M3M4) = {UlU 2 ){UsU 4 ) + {UiU'i) {u 2 Ui) + (M1M4) (M2M3) • 


(A3) 


So if these second order correlation terms are replaced by energy spectrum expressions as 
follows: 


A'Kk‘^{vi{k)vq{k')) = Pig{k)Ev{k)5k+k'fl, 

A'Kk‘^{Bi{k)Bg{k')) = Pig{k)EMik)6k+k'fl, 

A7re{B,{k)vg{k')) = P,g{k)HBv{k)6k+k'fl. 
kk 

(p. = A. _ .^j 


(A4) 
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Eq. 


can be rewritten like 


d 


+P^ + r'^)] {vq{p)v„,{r)vi{-k))= 


2Mins{-k) y^(47rp^) ^(47rr^) ^ Pnq{p)Pms{r)[Ev{p)Ev{r) - HBv{p)HBvir)] + 

p,r 

2Mqns{p) 5^(47rA:2)-'(47rr2)-ip,„(A;)P^,(r)[Ey(A:)Ev(r) - HBv{k)HBv{r)] + 

p,r 

2M^„,(r) 5^(47rfc2)-i(47rr2)-ip,„(fc)P,,(p)[Pv^(fc)Py(p) - HBv{k)HBv{p)] 


p^r 


= Liqmik,P,r-,t). (A5) 

However, when the fonrth-order correlation is decomposed into the combinations of second 
order terms, the snmmation of decomposed ones, i.e., right hand side of Eq. flASp . can be 
larger t han its actna l val ue. This can cause a negative energy spectrum, which cannot be 
allowed f Oguralll963 h So (Orszag 1970 1 introduced a n eddy damping coeffi cient of which 
dimension is l/t\ Its more detailed expression fjPonqnet et al.l 119761 ) can be contrived, 
but the dimension 1/t’ does not change. We assume it to be sort of a reciprocal of time 
constant for simplicity in this paper. Then, with a simple integration we can End the third 
correlation term: 


We can also calculate the representation of ‘A3’ in the same way: 


(A6) 


(A7) 


Since ‘A2’ and ‘A4’ are ‘-A1’ and ‘-A3’ respectively, Eq. flTTl) are 


{^ + 2iyk‘^){vi{k)vi{-k)) = ^ 2Miqm{k) 

p+r=fc 




(A8) 


It or L'’"^{k,p, r;T) is larger than (v(k^ + p'^ + r^) + ' or (vk^ + 

this equation can be markovianized. Thus, with the dehnition of a triad relaxation time 0(t) 
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flFrisch et al.lll975l) we get 
. d 


'dt 


2vk‘^)Ev{k) = 


A'Kk^Miqm{k) 

p+r=k 

_ g iyk +/l'fcpr)^ 

uk ~\- fJ^kpr 


]_ _ +P +?’ )+/^fcpr)i 

Z/(fc2 +p2 + ^2) 


L1qL{k,P,r]t) 


L'lqm{k,P,r;t) 


p-\-r=k *- 


(A9) 


Using a trigonometric relation: k ■ p = kpz, k • r = kry, p ■ r = —prx (Fig,5(a)) and 
dp dT=^^^dp dr, we can simplify this expression: 


1 

2 


l3 

dpdre''^"^(t)—{l - 2y'‘‘z‘ - xyz)\Ev(p)Ev{r) - HBv(p)HBv(r)]- 


(AlO) 


Eq. flT^ . flT^ can be derived using a similar way. 
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& Ky (y = 0.015. b^^^=0.0001) 




0 20 40 60 80 100 

time(X0.1, 288^) 


(a) 


(b) 


Fig. 1.— (a) Normalized Eyit) and EMif) for the incompressible flnid {Ptm —t cxo, resolntion 
256^ and 512^). (b) Normalized Ey{t) and Ej^it) for the compressible flnid {Ptm = 75 and 
Ptm = 7500, resolntion 288^). Time scale is contracted by 10% (xO.l). 



(a) 


(b) 


Fig. 2.— Energy spectra of incompressible flnid {Ptm —t cxo) (a) Ev{k) and EM{k) in the 
early time regime {t < 10.5) (b) Eyik) and pMik) in 24.9 <t< 32.6. As Em grows, large 
scale kinetic energy is transferred to magnetic eddies, bnt kinetic eddies in snbviscous scale 
receive energy from the magnetic eddies. 
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k (Pi'm^oo, 256^) 

(a) 




k (Pfh^oo, 

(b) 


512®) 
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(c) (d) 

Fig. 3.— (a), (b) Saturated Ey{k) and EMik) with resolution 256^ and 512^ for the in¬ 
compressible fluid, (c), (d) Saturated Ev{k) and Em{E) for the compressible fluids with 
resolution 288^. Saturated energy level of Ptm = 75 is higher than that of Ptm = 7500. 
Bottle neck effect appears. 
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Ev>Kk^° (y = 0.015) 



10 “ 


10 “ 


10 “ 


10 “ 


10 “ 



(a) (b) 

Fig. 4.— (a) Compensated energy spectrum k’^Eyik). (b) The flat reference line 

k^ means in EM{k), and slanted line k~^'^ is the scaling invariant line k~^. 


Em/Ev (u = 0.015) 



(a) 

Fig. 5.— The ratio Em/ Ey of hyper diffusivity increases with resolution. In case of physical 
diffusivity the equipartition of Em and Ey appears. 
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Fig. 6.— (a) The summation of three wave numbers should satisfy the condition p+r+k = 0 
(b) Upper triangle is the case of large p (small r): /c ~ p 3> r, a; ~ y, 2 : ~ 1, and lower one 
is for large r (small p); /c ~ r 3> p, a; ~ 2 ;, p ~ 1. 





